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Small-Signal  Gain  of  a  Free  Electron  Laser  in  a  Resonator  Gaussian  Mode 

Ilario  Boscolo  and  Juan  Gallardo 
Quantum  Institute 

UCSB,  Santa  Barbara,  California  93106 

ABSTRACT 

We  present  an  analytical  expression  for  the  small-signal  gain  of 
a  Free  Electron  Laser  (FEL)  in  the  presence  of  a  gaussian  mode.  To 
describe  the  electron  beam  evolution  we  use  the  one-dimensional  (1-d) 
Vlasov  equation.  Our  perturbative  result  is  valid  for  small  values  of 
the  parameter  q  (length  of  the  undulator  L  divided  by  the  Rayleigh 
range  . 
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I.  INTRODUCTION 


The  1-d  theory  of  FEL  has  been  thoroughly  studied^  and  by  now  we  have 
a  good  body  of  theoretical  knowledge  to  be  used  In  the  design  of 
actual  experiments.  All  these  papers  assumed  a  plane  wave  as  Input  field, 
as  differentiated  with  the  experimental  situation  where  we  have  gausslan 
modes.  This  Is  not  a  trivial  approximation  considering  that  the 
bunching  In  the  electron  beam,  and  consequently  the  loss  of  gain  of 
energy  by  It,  depends  critically  on  the  relative  phase  between  the 

electrons  and  the  electromagnetic  wave. 
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Colson  and  Elleaume  have  published  a  complete  treatment  of  this 
problem  for  several  magnet  designs  Including  the  varying  phase  and 

3 

amplitude  characteristics  of  a  gausslan  mode.  Also,  recently  analytical 
expressions  of  the  gain  have  been  obtained  using  the  single  particle  theory 
as  described  by  the  Lorentz  equations. 

4 

In  this  paper  we  used  the  Vlasov  equation  to  find  the  evolution 
of  the  electron  beam  In  the  combined  fields  of  a  copropagating  azimuthal 
symmetric  gausslan  mode  and  the  undulator  magnet.  The  radiation  field 
Is  described  by  the  Maxwell  equations  In  the  slow  amplitude  and  phase 
approximation^  (S.V.P.A.). 

The  FEL  dynamics  Is  properly  represented  by  the  coupled  system  of 
self-consistent  Vlasov  and  Maxwell  equations.  We  solve  It  using  pertur¬ 
bation  theory  In  powers  of  the  amplitude  of  the  Input  field. 


II .  THEORY 


The  single  particle  equations  of  motion  of  an  electron  in  the  combined 
wiggler  and  co-propagating  gausslan  radiation  field  are 
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so  the  permanent  magnet  field  ^  determines  Che  electron  trajectory  and 
the  phase  of  the  gausslan  field  the  energy  exchange. 

The  magnetic  field  of  Che  permanent  magnet  is 


B  ■  e»B  sin  k  z 
2  o  o 


Che  fundamental  gausslan  field  is 
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is  Che  beam  waist,  u)^  is  Che  spot  size  at  half-way  Che  wiggler;  and 
is  Che  Rayleigh  range;  the  phase  $  Is 
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is  the  curvature  radius, 
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Proceeding  a  la  Colson  we  introduce  the 
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the  dimensionless  time 


ct 


(0  <  T  <  1) 


and  the  velocity 


dc 
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We  also  define  the  total  phase  variable  as 

If  •  c  +  if(r,2) 

In  terms  of  these  quantities  the  single  particle  equations  of 
motion  read 

If  *  a(z)  cos  if  +  if 

•  • 

U  =  if  ■  V  +  (f 

where  u  is  the  generalized  velocity  and 
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Is  the  dimensionless  field  amplitude.  To  these  we  append  the  Vlasov 
equation  expressed  In  terms  of  t,  u.  and  ii  variables 
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Nov  we  turn  to  the  Maxwell's  wave  equation,  that  governs  the 
evolution  of  the  radiation  field  in  the  presence  of  the  electron  current 
described  by  Eqs.  3a)  and  5b). 

Employing  the  S.V. A. F. -approximation  and  following  Colson's 
procedure,^ the  resulting  wave  equation  can  be  written  as  a  pair 
of  equations  for  the  amplitude  a(r,z,t)  and  the  phase  ip(r,z,t), 
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The  dimensionless  current  density  Is 
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where  n^  the  particle  density  In  the  electron  beam.  The  complete 
set  of  equations  chat  describe  self-consistenCly  Che  physics  of  FELs 
Is  given  by  Eq.  5  and  6.  At  this  point  ve  restrict  ourselves  to  small 
signal,  low  gain  regime  where  the  amplitude  and  phase  of  the  optical 
waves  seen  by  the  electrons  do  not  change  with  time  except  for  the 


"Intrinsic  variation",  due  to  the  gaussian  nature  of  the  radiation 
field.  We  wish  to  point  out  that  it  is  here  where  we  depart  from 
the  plane  wave  calculation.  Therefore, 


where  we  have  neglected  the  transverse  r-dependence  of  the  gaussian 
field. 

III.  PERTURBATIVE  SOLUTION  OF  THE  VLASOV  EQUATION 

We  seek  a  perturbative  solution  of  Eq.  3b  in  powers  of  the 
amplitude  of  the  input  optical  wave  a(z).  In  first  order  we  have 

p(i/),u,t)  »  p^(p)  a(T)p^^^  (8) 

Since  the  coefficient  of  the  expansion  is  not  constant,  this  problem 
should  be  more  properly  treated  as  a  time-dependent  perturbation 
problem.  However,  we  observe  that  in  the  small  signal,  low-gain 
regime,  a(T)  can  be  assumed  a(T)  <  1  and  to  be  a  slowly  varying  function 
during  the  interaction;  hence,  we  are  justified  to  take  the  average  of 
the  amplitude  over  the  Interaction  length  <a(z)>  as  the  small  parameter 
in  the  perturbative  expansion.  Then  Eq.  8  can  be  rewritten  as 

P(<I»,U,t)  -  P_(U)  +  <a(T)>  p^^^  (iji.U.T) 
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Substituting  into  (5b)  we  find  for  the  perturbative  term,  the  equation 
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The  solution  is 
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Using  the  second  relation  in  Eq.  (5a),  we  can  perform  the  integral 
obtaining 
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The  electron  distribution  function  in  first  order  in  <a(z)>. 


IV.  SMALL  SIGNAL  GAIN 

By  definition  the  small-signal  gain  is  the  relative  power  variation 


/  a^  dS  -  /  a?  dS 

■'  o _ ■'  1 

f  2 

/  dS 


(11) 


where  a^  and  a^  are  respectively  the  dimensionless  output  and  input  field 
To  calculate  the  final  field  a^  we  use  Eq.  (6)  obtaining 

1 

a  ■  a.  -  j  [-  dx  <cos  ')>>  .  (12) 


Ue  assume  Che  initial  electron  beam  to  be  monochromatic  in  energy;  therefore 


the  zero-order  distribution  function  must  be  a  constant,  sharp-defined 
function  of  Che  generalized  velocity  y  along  the  interaction  region. 

We  write  the  initial  electron  distribution  as. 
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with 
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where  we  have  approximated  caking  the  average  value  of  Che  time 
derivative  of  the  gausslan  phase  <<J>  over  the  length  of  the  wiggler. 

This  approximation  is  consistent  with  our  previous  discussion  for 
handling  Che  time-dependent  gausslan  amplitude.  Since  we  know  the  electron 
distribution  function,  it  is  straightforward  to  compute  the  phase 
distribution 


<COS  l|)> 


dip  dy  COS  4)  p(4),u,t) 


(15) 


The  variation  in  Che  field  amplitude  due  to  Che  perturbed  part  of  the 
electron  distribution,  is  6a;  hence  we  can  write 
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After  performing  the  integration  we  get 
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where  g(y^)  is  the  usual  lineshape  function  for  the  gain.  We  recall 
that  the  variation  of  the  field  6a  is  defined  over  the  transverse  area 

of  the  electron  beam.  Combining  (16)  irtth  (11) ,  we  find  for  the  optical 
gain 


2a.  5a  v 
—i - LI 

f  dS  a? 


Introducing  explicitly  the  electron  current  I 
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the  averaged  amplitude 
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and  the  input  flux  *•  1/q.  We  obtain  the  final  expression  for  the 

small  signal  gain. 
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We  point  out  the  following  features  in  Eq.  20, 

a)  the  functional  form  of  the  gain  is  given  by  8(Uq) »  the  usual  lineshape 
function;  the  gaussian  beam  only  introduces  a  shift  of  the  entire  curve 
towards  higher  values  of  the  detuning  parameter  \)^,  as  shown  by  Eq.  14. 

b)  the  function  in  square  brackets  is  a  slowly  varying  function  of  q 
(Figure  2)  reaching  a  maximum  for  q  -  3.  This  value  of  q  will  optimize 
the  design  of  an  FEL  resonator.  The  shift  of  the  gain  curve  toward 
higher  values  of  is  approximately  linear  with  q. 


We  remark  that  the  classical  gain  formula 


in  Eq.  20  is  the  result  of  our  approximation  of  replacing  the  amplitude 
and  phase  of  the  radiation  field  by  their  averages  <a(z)>  and  <0(z)> 
over  the  length  of  the  undulator. 

These  results  coincide  with  those  of  references  2  and  3.  In  addition 
in  those  papers  it  was  shown  chat  the  gain  curve  is  slightly  distorted 
and  that  the  absorption  maximum  is  larger  than  the  gain  maximum. 


V.  CONCLUSIONS 


We  have  presented  an  analytical  expression  for  the  FEL  small-signal 
gain  in  the  presence  of  a  gaussian  mode.  The  gain  curve  is  shown  to  be 


shifted  and  distorted  with  respect  to  the,  by  now,  classical  1-d  formula. 

The  maximum  gain  Increases  (linearly  for  q  0)  with  q  reaches  a  mavimniti 
and  then  decreases;  q  ~  3  Is  the  optimum  value  to  be  used  In  the 
design  of  FEL  resonators.  The  shift  of  the  gain  curve  is  also  linear  with 
q  and  changes  substantially  Its  rate  of  change  at  q  ==  3.  All  these  results 
are  In  fairly  good  agreement  with  previous  works. 
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